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The diffractive f-slope, Bv, of elastic charmonia (V =J/if), ip') photo- and leptoproductions off a 
nucleon is studied at low |t| (< f GeV 2 ) in the leading logarithmic approximation of perturbative 
QCD, with a special emphasis on the space-time evolution of the cc-dipole. We obey a framework 
based on QCD factorization, which describes a certain Fermi motion effect due to the c(c)-quarks 
in the proper manner P, Q and includes appropriately kinematical corrections of the momentum 
transfer Ax in the t-channel. Assuming the universal two-gluon form factor of the nucleon, we 
show that the difference of f-slopes for J/ip and if is dominated by the contribution from the 
dipole-charmonium transition process. The calculated difference is found to be Bj^ — B^i ~ 0.53 
GeV -2 for the photoproduction, in agreement with HERA data. We also calculate the i-dependence 
of the total cross sections for several center-of-mass energies W . A good agreement of the results 
with the available data demands that the mass scale appearing in the gluon form factor should 
significantly decrease with increasing W. 
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I. INTRODUCTION 

Recent measurements of high-energy diffractive photo- 
and electroproductions of charmonia off a nucleon, 7^*' + 
N(P) -> V + N(P') (V = J /if, if'), are unique sources 
of information about both structures of the target nu- 
cleon and the production process of charmonia. The 
t-distribution of differential cross sections, where t is a 
square of momentum transfer to the vector meson de- 
fined as A = P — P' , i.e., t = A 2 ~ — Aj_, provides us 
with information concerning the spatial structure in the 
plane perpendicular to the photon- nucleon reaction axis. 
In this case, one can study the impact parameter distri- 
bution of the gluon in the nucleon, as well as the space- 
time evolution of the dipole-like cc state, which is cre- 
ated initially by the photon's fluctuation and hadronizes 
to charmonia after a direct scattering with the gluons in- 
side the nucleon, in the transverse space. The former is 
never accessible in the well-known inclusive deep- inelastic 
scattering, from which one can extract only parton dis- 
tributions of the nucleon with respect to the longitudinal 
momentum fraction. In the limit of large photon virtu- 
ally Q 2 , a remarkable prediction by QCD factorization 
tells us that the t-dependence of a production cross sec- 
tion is solely determined by a universal two-gluon form 
factor of the nucleon, independent of the type of vector 
mesons produced 0. 

On the other hand, the effect of the dipole dynamics 
on the ^-distribution of cross sections is considered less 
important than the one due to the nucleon structure, be- 
cause the cc state is highly squeezed at large Q 2 and, 
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naively speaking, can be regarded as an almost point- 
like state during the whole process. In the diffractive 
charmoniaphotoproductions observed at the HERA ex- 
periment (J, lil 111 j however, only a hard scale is provided 
by the charm quark mass m c , e.g., Q 2 ~ to 2 . In this case, 
the initial cc state is dominated by transverse polariza- 
tions of the quasi-real photon and thus should not be so 
squeezed. Therefore, the evolution of the cc state in the 
transverse space most likely gives a non-negligible contri- 
bution to the ^distribution of the production amplitudes. 
It is a significant work to estimate how the ^-distribution 
of cross sections is sensitive to the structure of the cc 
state as functions of Q 2 and W. 

For this purpose, it is natural that we formulate the 
cross section at finite t along the familiar frameworks 
based on QCD factorization Q, 0, @ , which give a good 
description for the diffractive J/tp photo- and electro- 
productions. The perturbative QCD analyses (pQCD) 
for the J I if) photoproduction at t = suggest a strong 
suppression of the total cross sections due to the Fermi 
motion of the c(c)-quarks, especially relative transverse 
motion between c and c. In particular, it is more pro- 
nounced for a radially excited if'(2S) production, because 
the radial size of the cc-dipole in the final hadronic state 
is about twice as large as that of J/ij) and then the trans- 
verse motion of the c(c)-quarks could be more active ow- 
ing to the existence of a node in the if' wave function 
lj. Thus, detailed information concerning the internal 
structure (i.e., motion or radial size) of the cc state in the 
dipole-charmonium transition process are considered to 
be inherent in the cross sections of the diffractive J /if and 
if' productions. One possible method to extract such in- 
formation suitably from these processes would be to cal- 
culate the i-distributions of the cross sections. This fact, 
further, motivates us to find out information associated 
with only the dipole structure of the cc state by making 
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use of both the J/ip and ip' production cross sections. 
This would be carried out in the following way: experi- 
ments at HERA have observed the i-distribution inten- 
sively for the J / ip photoproductions 0, and recently, 
also for the ip' case 0. Each i-slope, By, extracted as 
a result of an exponential fit to the data, exp(i?yt), in- 
cludes both information on the nucleon and the cc state, 
as mentioned above. We assume that the t-slope is of a 
simple form described by a sum of a universal nucleon 
form factor and the contribution from the dipole state at 
large W (see Sec. II A). Therefore, a difference of i-slopes 
between J / ip and ip' is free from the contribution of the 
universal nucleon structure, and is determined only by 
the finite size effect of the cc-dipolc. This difference is in 
fact calculated as the convolution of the dipole scattering 
amplitudes with the wave functions of J/ip and ip' . 

Motivated by these interests, in this paper, we calcu- 
late a differential cross section as a function of \t\ (< 1 
GeV -2 ) in diffractive (elastic) charmonia (J '/ip and ip') 
photo- or electroproductions off a nucleon. We employ fa- 
miliar QCD factorization formulae with the helicity rep- 
resentation in the Brodsky-Lepage approach for hard ex- 
clusive processes 0. It is calculated under the leading 
logarithmic approximation (LLA) of pQCD, which is rea- 
sonable for the interaction of a small transverse-size cc- 
dipole. A reliable prediction for the i-slopes requires a 
sophisticated pQCD analysis for the space-time evolution 
of the cc-dipole, particularly giving a precise description 
of the Fermi motion between c and c. On thispoint, our 
approach in the dipole picture developed in 0, could 
be more appropriate than those in Ref. [H], because 
our formulation deals with 0(v 2 ) corrections due to the 
Fermi motion effect in a proper manner, where v is an 
average velocity of the (anti-)charm quarks in the char- 
monium rest frame. The previous approaches Il0j | are 
not satisfactory in view of reliable estimation of the Fermi 
motion effects, which should take into account a proper 
projector onto the S-wave wave function , the con- 
tribution of off-shcllncss in the spinor matrix elements of 
cc state and other 0(v 2 ) corrections 0. The approach 
of is especially inadequate for the description of ip', 
because their light-cone wave functions (LCWF) of char- 
monia are constructed in an oversimplified fashion, such 
that the dependence of the ratio of ip' to J/ip cross sec- 
tions on Q 2 at t = fairly underestimates the HERA 
data at small Q 2 [1 [U El . 

Another development of our model is to include a kine- 
matical correction from finite Aj^ due to a recoil of the 
target nucleon, which should explicitly appear in the ar- 
guments of nonperturbative charmonium wave functions 
(a, fcj_, A^), with a and k± being the longitudinal 
momentum fraction and the transverse momentum of the 
c-quark, respectively. Previous works [3, 1 1 0| simply take 
the charmonium wave function \E , y(a,fc_L) in the limit 
of Al = and miss kincmatical corrections arising from 
the finite Aj_ . Since the fc^-dependence of the wave func- 
tion is dominated by the region, where |fcj_| is less than 
the inverse of typical hadronic size, ~ 0.5 GeV, the con- 



tribution of |A^| is not small compared to |fej_|, even if 
we restrict the region of \t\ = A^_ to less than 1 GeV 2 
in our analysis. We explicitly deal with such a correction 
under the LLA. 

With these improvements, we derive the i-slopes of 
both J / ip and ip' productions and investigate the Q 2 and 
W dependences respectively. We find that the contribu- 
tion of the cc state to the i-slope is not negligible and 
the sizes reach about 5 ~ 10 % of the whole, which be- 
come largest at Q 2 — for the J/ip production, and at 
Q 2 ~ 10 GeV 2 for the ip' one. The difference between 
t-slopes of J / ip and ip' is then calculated to study the dy- 
namics of the cc-dipole. We emphasize that this quantity 
provides us with an opportunity for direct comparison of 
the pQCD calculation of the dipole contribution with the 
experimental data. The result is actually consistent with 
the data from HERA at Q 2 = 0H|(f. 

We also calculate the i-dependent differential cross sec- 
tion of J/ip by combining the dipole part and the univer- 
sal gluon form factor of the nucleon. We fix a value of 
the mass scale appearing in the dipole form of the gluon 
form factor so as to reproduce the observed i-dependence 
0| . Precise determination of the mass scale should be 
possible, because the present work already gives a real- 
istic description of the dipole contribution to the cross 
section. In order to obtain a good fit to the data of In- 
dependence, it turns out that the mass scale of the gluon 
form factor should become significantly smaller with in- 
creasing W . 

The paper is organized as follows: in Sec. II, we give a 
detailed description of our formulae in the LLA, focusing 
on the Aj_ dependence. Sec. Ill contains our numerical 
results of the t-distribution coming from the cc-dipole 
in the J/ip and ip' productions. Here, the Q 2 and W 
dependences of each i-distribution and their difference 
are shown and compared with HERA data at Q 2 = 0. 
Also, the VF-dependence of the mass scale of the form 
factor is discussed through fits to data. In Sec. IV, a 
summary and discussion are presented. 

II. FORMULATION 

A. Description of t-slope in QCD factorization 

In the pQCD analysis, the corresponding diffractive 
amplitude at high W can be interpreted as a sequence of 
several steps separated in time, as demonstrated in . 
It schematically has a factorized form 

M = $>y ® Accgg ® f ® * 7 . (1) 

Here <3/ 7 is the LCWF of a photon describing the cc fluc- 
tuation from the photon, A C cgg is the hard scattering 
amplitude of the cc-pair off the nucleon via two-gluon 
exchange in the i-channel, T represents the t-dependent 
gluon distribution in the nucleon and tyy is the char- 
monium LCWF including the soft hadronization pro- 
cess of the cc state. The calculation for the transition 
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-y(*) V is based on a technique developed in 0, 0, 
which accounts for moderate sub-leading effects to 0(v 2 ) 
due to the Fermi motion of the cc state, especially impor- 
tant for ip'. On the other hand, the process-independent 
gluon distribution T(xi,X2, l±, t) in the nucleon is as- 
sumed to be a product of the off-forward (skewed) unin- 
tegrated gluon distribution f(x±, X2, l±) with a universal 
two-gluon form factor T(t), i.e., T = f(x\,X2,l\) X T(t) 
[l3L III ]. Here, x\ and X2 denote the longitudinal mo- 
mentum fractions of outgoing and incoming gluons re- 
spectively, and l± the transverse momentum of the out- 
going gluon (see Fig. QJ. The form factor T(t) has the 
same expression as a usual electromagnetic one of the nu- 
cleon, so called "dipole form" T(t) = 1/(1- t/M 2 ) 2 with 
the mass scale M 15]. Here, it should be noted that 
the cross sections of these processes are written as the 
quadratic form of the nucleon form factor multiplied by 
the charmonia LCWF. Thus, the diffractive charmonium 
productions could be sensitive to those nonperturbative 
objects. 

The experimental differential cross section da/dt cx 
\M\ 2 is usually parametrized as exp(Bvt) with the 
diffractive slope By, as already mentioned. The t-slope 
is extracted by the logarithmic derivative of the cross 
section over t. Owing to the factorization ansatz, the 
calculated i-slope can be completely separated into two 
contributions of the nucleon part, B^lt), related to the 
nucleon form factor T(t), and the dipole part, By P (t), 
associated with the transverse evolution of the cc-dipole 
in the transition process 7^*' — > V: 

B v (t)=B N (t)+B^ p (t), (2) 

where B N (t) = 4/(|f| + M 2 ). 

As a result, the difference of i-slopes between J/ifj and 
ip' productions is independent of the universal nucleon 
form factor and thus equals the difference of the dipole 
parts, 

(3) 



Bj^{t)-B^{t) = B a ;f^t)-B^{t) 



This expression enables us to investigate the structure of 
the cc dipole in detail through a direct comparison with 
the data. The naive geometrical interpretation may sug- 



gest that the difference £> 1S P — B 2 g is negative, because 
the smaller the size of color-singlet objects (or qq pair), 
which is approximately the inverse of mass of the vector 
meson, the smaller the ^-slopes [Tij . As emphasized for 
charmonia productions in |l(J, however, when one con- 
siders the production of the radially excited state such 
as ip'(2S), such a naive expectation is broken because of 
the node effect in the 2S wave function. Actually, recent 
HERA data at the photoproductions seem to support 
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N(P) -» V(q + A) + N{P — A) (V = J/ip,ip'), where 
P and q denote the four momenta of the photon and 
the incoming nucleon, respectively, and A the momen- 
tum transfer in the t-channel, as illustrated in Fig. ^ 
The tot al ccnter- of-mass energy of the j^-N system, 
W = WjP + q) 2 , is then assumed to be much larger 
than the photon's virtuality Q 2 and the heavy-quark 
mass m c . Following the usual fashion, we perform the 
standard Sudakov decomposition of all momenta labeled 
in Fig. ^ by introducing two null vectors q' and p', i.e., 
Q = q'-(Q 2 /s)p>, P = (M 2 N /s)q'+p\ k = aq' + f3p' + k ± , 
I = nq' + xip' + l± and A = (t/s)q r + Sx p' + A±. Here, 
s = 2q' ■ p' ~ W 2 — M^r + Q 2 with the masses of the nu- 
cleon and charmonium, Mn and respectively, and 
= —t. We assume that q' + and p'~ are much larger 
than q'~ and p' + in our reference frame. The coeffi- 
cients 13 and n are of the order 1/s. x\ has the relation 
x\ = (M 2 S + Q 2 ) I s + 0(l±) with the squared mass of the 
intermediate cc state, M 2 5 = (k ± + m 2 c )/[a{l — a)]. The 
second term of x±, 0(l±), corresponds to the contribu- 
tion beyond the LLA, and thus we neglect this term here- 
after. The longitudinal momentum fraction of the incom- 
ing gluon, X2 — x\—8x, is given as X2 — (M 2 E —M 2 +t) / s. 
Since X2 « 11 in the small region of \t\ < 1 GeV 2 , this 
means the off-forward (or skewed) kinematics for the two- 
gluon exchange. 




FIG. 1: Four Feynman diagrams for 7'*' + N — ► V{= 
J/il),i[)') + N with momentum labeling. 



The differential cross section with the small momen- 
tum transfer t to the nucleon has the form 



da(^*~>N -> VN) 



dt 
1 
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J2 \lmM poL (t) + i ReM poL (t) \ 



pol.=L,T 



B. Differential cross section in the LLA of pQCD 



(4) 



In the LLA, we formulate the diffractive photo- or elec- 
troproduction of charmonia off the nucleon, j^*\q) + 



where the suffix 'pol.' denotes the polarizations of the 
incoming photon. The imaginary part of the amplitude 
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M poL is of the form 
lmM poL {t) 




dl 



± ll(l±-A x ) 2 



/2N C Jo J 16?r 3 
x T( Xl ,X2,ll,t)IP° L ^ v \a,k ± ,l ± ,A ± ) (5) 

with the number of color N c = 3 and a s = g 2 /(4ir), 
where jp oL (i^ v ) represents the transition 7W — > V such 
as 



I pol -^ v \a,k ± ,l ± ,A ± ) 
= \% A H<xM-% B H<*,kx,lj 



XX' 



-^(a,k ± ,l x ,A A 



A' A 



X * 



VAA 



(a,k ± , A ± ) 



(6) 



with the helicities of charm and anti-charm quarks, A' 
and A, respectively. The function T(xi,X2,l±,i) describ- 
ing the lower part of Fig. ^ as mentioned in Sec. II A, 
has the form T = f(xi,X2,l 2 ±)T(t) with the two-gluon 
form factor, 



r(t) 



(i-t/O' 



(7) 



where 77i2 g denotes a mass scale relevant to this pro- 
cess 0. The skewed unintegrated gluon distribution 
f(x\,X2,i 2 ±) is related with 



f(x 1 ,x 2 ,l 2 l ) = 



dx 2 G(xi,X2,l 2 L ) 
<91n I 2 , : 



(8) 



using the skewed gluon distribution G(xi, x%, l ± ). The 
effect of skewedness can be then incorporated into a con- 
stant enhancement factor R at small X\ like the HERA 
experiment |16(. Therefore, we rewrite the function 
f(x\,X2,l±) in terms of usual diagonal gluon distribu- 
tion as 



f(xi,x 2 ,l ± ) = 



2 , Rd Xl G( Xl ,l\) 



din I 2 



(9) 



Here, R = 1.1 ~ 1.2 for the J/ip photoproduction 0] 
and this factor gives an enhancement by 20 ~ 40 % for 
the cross section. Thus, we express the imaginary part 
of amplitude in terms of the diagonal gluon distribution, 
where hereafter we denote x\ as x for simplicity. Then, 
the real part in Eq. (jJJ is related to the imaginary part 
as ReA'l = — (7r/2)(9ImA / J)/(91nx) in the perturbative 
analysis 0- 



C. Light-cone wave function of photon 

First, let us construct the LCWF of the incoming pho- 
ton ^ 7 in Eq. 10, based on the LC perturbation theory 



00 It is defined as the following spinor matrix element 



7A'A 



UX' (gt 

Ja. 



jpol. 



vx(qi) 



-ee c a(l 



a/I - a k ( f> 2 + ml + a (l - a)Q 2 ' 



(10) 

where e poL denotes = (q' + xp')/Q for the longitudinal 
polarization, ande^r = (0,0, e^) withe]^ = t^(1i7*) an d 
7 = ±1 for the transverse polarization. Here, • q = 
e!^ • — 0. e c is the charge of charm quark 2/3 in the 
units of e. and are the momenta of charm and anti- 
charm quarks with the helicities A' and A respectively, 
for each diagram (i) illustrated in Fig. ^ We define fej^ 



fe 



■fx, fe 



(C) 



fc I -I I 



■ A 1 and 



k^ = k± + Aj_ for the corresponding diagram. After 
calculation of the spinor matrix elements, the results of 
the longitudinal polarization are given by 



long. 

7 A'A 



f: 



for A' = A, and 



long. 

7 A'A — 



2a(l - a)Q 2 



k { l )2 + m 2 c + a{l - a)Q 2 



(11) 



(12) 



for A' = —A. Here, the factor 1 in the first term of Ea. l|12|) 
disappears, when one takes a sum of four diagrams in 
calculation of the amplitude, and thus one neglects the 
factor. The results of the transverse polarization are 



T trans. ( 7 ) 
^ 7 A'A 

for A' = A, and 

- trans. (7) 



\/2ee c XS-, 



xm c 



.«2 



+ m 2 c + a(l - a)Q 2 



7 A'A 



for A' = -A. 



2ee e (q - <5A, 7 )el ■ fe 
'e ( l )2 + m 2 + a{\ - a)Q 2 



(13) 



(14) 



D. Light-cone wave function of charmonium 

Next, we describe the LCWF of the outgoing charmo- 
nium produced through the soft hadronization process 
of the cc state, "fy. It is defined in analogy with that 
of a photon. However, one should pay attention to new 
contribution from non-zero momentum of Aj_, different 
from the case of the photon. Namely, the center of mass 
motion of the charmonium deviates by A j_ to a direction 
perpendicular to the incoming photon, due to a recoil of 
the target nucleon. Therefore, we introduce new null 
vectors q' and p' to redefine (q + A) + as the longitudi- 
nal direction of fy Here, q' and p' are expressed as 
q' = q' - (t/s)p' + A_l and p' = ~(t/s)q' + p' - A± up 
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to 0(1/ s) in terms of q' and p'. The momenta of two 
quarks in the cc state, k + A and k — q, are written as 



k + A = ag' 



at , Q 2 -i + M 2 



+fex + (1 - a)Aj_, 
fc — q = (a — l)q' + 

+k± + (1 - a)A 



at 
s 



- t 



(15) 



Using these relations, we find that the relative momen- 
tum between c and c is fcx + (1 — a)Aj_. The polarization 
vectors of the charmonium are also rewritten in terms of 
q' andp' as = [q'- (M 2 +t)p'/s + A±]/M^ for the lon- 
gitudinal polarization, and £+ = e^+2(A±-e^)(q' — p')/s 
for the transverse one. Here, e° • (q + A) = ■ (q + A) = 
e^p- = 0. Using those polarization vectors, the LCWF 
of the charmonium $y has the form 



*£ft(a,fcx,A x ) = 



<^(a,fcx + (l-a)Ax 



A) 



(16) 



where 7?. represents the projector introduced in P, 
which keeps the cc state to be low-lying IS state, given 
as K = [1 + (4+ A)/M^]/2. Following 0, we further 
take into account the off-shellness of the spinor matrix 
element, v\^.° *1Zu\i , as the proper treatment. This 
spinor matrix element representing the effective "cc — * V 
vertex" may be off the energy shell and thus the total 
energy carried by the two spinors v\ and uy is generally 
different from the energy of the vector meson. This con- 
tribution should give 0(v 2 ) corrections in the final result 
(In more detail, see @). For the longitudinal polariza- 
tion, the results lead to 



long.* 



vw = -V2A(1 - 2a)e [ l ] ■ (fej 
(a, fcx + (1 - a) Ax 



(l-a)Aj 



2a(l - a)M^ 



(17) 



for A' = A, and 



long.- 
VW 



(fcx + (1 - a) A ± ) 2 + m c (M^ + m c ) 



a(l - a)M 2 



for A' = A, and 



trans. (7 )* 
VW 



2el*-(fex + (l-a)Aj 



x [(a-5 7!A )(M^ + 2TO c ) + (l-2a)m c ] 
0;(a,fcx + (l-a)Ax) 



2a(l - a)M 2 



(20) 



for A' = —A. The suffix 7' denotes the polarization 7' = 
±1 in the transverse direction. 



E. Detailed derivation of differential cross section 

We focus on the process of S'-channel helicity conser- 
vation (SCHC) between the initial photon and the fi- 
nal charmonium, because the HERA experiment demon- 
strates that, to good accuracy, the hypothesis of SCHC 
holds in the kinematical range 30 < W < 200 GeV and 
\t\ < 1 GeV~ 2 [5j. We obtain the following expressions 
f or jpoi-h^v) of Eq ^ using Eqg ||nj-J2J and (fTVfi - 

(|20|l : for the longitudinal channel of j(L) — » ip(L), 

I L ^ L = 2ee c a(l - a)Q $1 

(fcx + (1 - a)A ± ) 2 + m c (M^ + m c ) 



1 



a(l - a)M? 



4, 



(21) 



x ^(o, fex + (1- a) A x ) ; 
for the transverse channel of j(T) — > ip(T), 

jT( 7 )->T( 7 ) = ee c 

a(l - a)M 2 
x [m c {(M v ,+m c )m c + a(l-a)M|}$i 
- 2 {(2a(l - a) - 1)M 4 , - m c } e\ • $ 2 



x el* ■ (fex + (1 - o)Ai)] 
x^(a,fci + (l-a)Ax) . 



(22) 



Here, conveniently we defined the sum of energy denom- 
inator as 



4>i 



1 



1 



k\ + Q 2 



(fe_ 



Q 2 



1 



x 4>*Ja, fex + (1- a) Ax) 



(18) 



for A' = —A, where = ^=(l,Ai). For the transverse 
polarization, similarly, 



$ 2 = 



(fex- 

fc 



—2 



(fex + Aj 



Q 



(23) 



fej 



k 



fel 



(fc ± + i ± ) 2 + Q 



(t> = y^A |(M^ + m c )m c + a(l - a)M 2 } 

+25y,_ A {el'* • (fex + (1 - a)A x )} 
4>;(a,fex + (l-a)Ax) 



(fc ± + A ± -Z ± ) 2 + Q 
^a(l-a)Q 2 . 



(fex + A ± ) 2 +^ 



(24) 



2a(l - a)M 2 



(19) 



where Q = mj 

Now we are in a position to perform analytical inte- 
grations over l± and k± in the amplitude (jSJ). First, we 
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explain the integration over l±, assuming the range of l 2 j_ 

relevant in this process as < l\ < Q +k\. This 
permits one to make an approximation — A^) 2 ~ 

l/l± for the gluon propagators in Eq. JSJ). Then, the 
Z_i_-integration for $1 in the integrand of Eq. (J3J) leads to 



f+fc " ,. a s (ll)dxG(x,l ± ) ^ 

dlj_ — ^ — $i(a, k±,l±, Aj 



dl 2 ± 



V 



Tra s (Q + k z L )xG{x,Q + k J ± ) F t {a, k ± , A ± ) 



(25) 



with 



Fi(a,k±,A±) = - 



2Q 



k\ + Q 2 



k 2 ± + Q 2 



2Q 2 



A,? 



Q 2 } {(fe_ 



A,) 2 



(26) 



Here, we integrated from the low energy cut-off /i 2 (^ 

A qcd) t0 Q 2 + k l HG1 over Similarly, for £ | ■ $ 2 
in the integrand of Eq. |J5J, we get 



a s (i 2 L ) dxG(x, I 



dl ± — T 



eX • $ 2 (a;,fc_L,«±, Aj 



7ra s (Q +fc 2 L )xG(a;,Q + fcl) eX ■ F 2 (a, fe±, A_l) 

(27) 



with 



el-F 2 (a,k ± ,A ± ) 



2Q 



el-(k ± + A ± ) 



Q' 



(28) 



where we neglected O ((£j_/ 'Q 2 ) 2 ^j terms in the LLA. 

Next, we analytically perform the angular integration 
over k± in Eq. J^J). We introduce a new integral variable, 
K± = k± + (1 — a)Ai, instead of k±, to rewrite the 
nonperturbative wave function (f>^(a,k± + (1 — a)Ai) 
as a function of one variable, K±. Here, we stress that, 
as mentioned in Sec. I, the contribution of A^ is not 
suppressed compared with that of k± in the argument 
of wave function, and thus we should exactly include the 
kinematical corrections of Aj_ by such a redefinition of 
the argument. When we define the angle 8 as K ± ■ A± = 
|-Kj_||Aj_| cos#, the ^-integrations of Fi and F 2 yield the 



following results: 



r- 2 TX 

\ d0F 1 (a,K ± ,A ± ) 
Jo 



4tv 



Kl+Q 2 



l , 2{(o 2 + (l-q) 2 )AX-Q } 
Kl+Q 2 
12(a 2 + (1 - a) 2 )AiQ 2 



K\ + Q 2 



12(a 2 + (l-a) 2 )AiQ 4 



(k\ + Q' 



(29) 



dO el-F 2 (a,K ± ,A ± ) el*-K ± 



AttKIq 2 
K 2 ± + Q' 



3(a 2 + (l-a) 2 )Ai 



Kl+Q 2 



6(a 2 + (l-a) 2 )AX 



Kl+Q 2 



(30) 



Here, we left the ^-dependent terms up to 

O (aI/(KI + Q 2 )}. Therefore, in the photopro- 

duction, this implies that the applicable range of t 
should be \t\ <C m 2 and we should discuss numerical 
results in the small regions of \t\ < 1 GeV 2 in Sec. III. 

Finally, using Eqs. (|29|l and l|3UI) . we arrive at the whole 
expressions for the imaginary part of amplitude 10: for 

7 (L)->V(£), 



lmM L (t) = - 



2V2w 2 ee c QsR f 1 



/Nr. 



daa(l — a) 



(Q +kl)xG(x,Q +kl) 
dK 2 ± <P*4a,K ± ) 



16tt 2 



(kI + q 2 )' 



( Kl +m c (M^ + m c ) \ 
{ "(l-a)M 2 ) 

2(a 2 + (l-a) 2 )(-t) ij ' 

Hi{a,Kx) H —5 H 2 {a,K 1 _) 

Kl + Q 



(31) 
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for 7 (T)->VCr), 

t ijT t .\ V2n 2 ee c sR 
Im M (t) = - 



da 



VKM 2 Jo a(l-a) 

(Q +k 2 ± )xG(x,Q +k 2 x ) 
dK 2 (f>*Ja,K x ) 



16tt 2 



Q- 



-m c {{M^ + m c )m c + a(l - a)M^} Hi(a,Kj_) 



+ 2 {(a 2 + (1 - a) 2 )M^, + m c } 
2(a 2 + (1 - af){-t) 



Kl + Q 2 



K 2 +Q 2 



ReM 



L(T) 



A 



R< 



+ (— t)A^f T) . Using these amplitudes 

2\ 



l L(T) 

to the first order of t, By P (W, Q 2 ) is given by 



B V * P (W,Q 2 



Aim Aim j_ A Re A Re 



= L,T 



(34) 



III. NUMERICAL RESULTS AND 
COMPARISON WITH DATA 

In Fig. |2 we show the results of the dipole contribu- 
tion to the slope parameters, By P , as a function of 

dip 

and solid line for their differ- 
85 GeV. B 



Q~. Those results are plotted as dashed line for B 
dash-dotted line for B2f ', 
ence B j P ,,—BJ p , at the fixed energy W 



J/4> 



x |-m c ((M^ + m c )m c + a(l - a)M|) ff 2 (a, K ± ) 

+3 ((a 2 + (1 - a) 2 )M^ + m c ) 
K\Q 2 



Kl + Q 



'■=sHi(a,K ± ) 



(32) 



Here, we define 



ffi (a, «■_!_) = 1- 



2Q 



H 2 (a,Ks_) = 1 



6Q 2 



6Q 



(33) 



To determine the nonperturbative part of the charmo- 
nium wave functions, 4>ip(a, K we use the wave func- 
tions obtained by solving Schrodinger equation with the 
realistic potential 0, Il2fl We adopt the Cornell poten- 
tial model with the corresponding quark masses m c = 1.5 
GeV |l7|. giving a good description of the charmonium 
wave functions. Here, we rewrite the non-relativistic 
wave function, <fi^ R , originally obtained in terms of three 
momenta, as a function of the LC variables, </>,/, (a, K 
by simple kinematical replacement [7lll2|. For the diag- 
onal gluon distribution G in Eqs. I|31l) . (|32[) . we employ 
Gliick-Reya-Vogt parametrization for the next-to-leading 
order fits [l8|. 

Using Eqs. I|31[) and H32|) . let us construct By P (t) in 
Eq. Pjl. which comes from the upper part describing 
the transition 7M — > V in Fig. One should note 

that, in our approximation up to the first order of t, 
B dlp (t) is a function of Q 2 and W, whereas indepen- 
dent of t. Conventionally, we denote the imaginary part 
of the amplitudes as ImM^T) — ^ut) + ( — *)^l(t)> 
where A 1 ^^ corresponds to the amplitude independent 
of t in Eqs. (|31|l and (|32|l . and then -Awm denotes the 
amplitude proportional to t. Similarly, for the real part, 
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FIG. 2: The slope parameters Bj,^, , uuu ^j/^ 
plotted as a function of Q 2 at the fixed energy W = 85 GeV. 
A blob at Q 2 = represents Bjm — B^i from the data @ for 
a reference. 



has the maximum value of B 



0.48 GeV" 2 and de- 
creases with increasing Q 2 , while Bj p has the negative 



value near Q 2 — 0, B 



- -0.07 GeV" 2 , and rapidly 
increases with increasing Q 2 , taking the positive value. 
This feature changing the sign of the i-slope with Q 2 
appears only in the tj)' production. This is an indica- 
tion of the node inherent in the radial direction of the 
2S wave function. At large Q 2 , however, they approach 
each other, and then their difference reaches almost zero, 
because the size of the cc-dipole state is squeezed with 
increasing Q 2 and the amplitudes are evaluated near the 
origin of the charmonium wave functions. We find the 
maximum value of the difference at Q 2 = 0, which is 



B 



dip 



B 



dtp 



_ "ty — 0-53 GeV 2 . Using Eq. J3J), we can di- 
rectly compare the result with the experimental data. In 
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fact, it turns out that our result is consistent with HI 
data 0, which is Bj/^ — B^ ~ 0.68 GeV -2 showed as a 
blob in Fig. [21 in the sign and the magnitude. 

Similar behavior of the i-slope difference as a function 
of Q 2 was also confirmed in the analysis of [lfj for the 
transverse polarization. Their results show that the dif- 
ference is about 0.25 GeV -2 at Q 2 = and W = 100 
GeV, which is about 2 times smaller than our result, and 
falls down more slowly with increasing Q 2 . 

In Fig. we plot the slope parameters B d j?, , B d f ', and 



B 



Bj p , as a function of W at Q 2 = 0, similar to 
Fig. [21 Apparently, both B dtp ^ and Bj p show an almost 
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3: The IF-dependence of the slope parameters Bjf^, 
and B rf.i. — B^f ', for photoproductions, compared with 



dip 

ZEUS data p. 



flat behavior as a function of W and their absolute values 
even decrease slightly with increasing W. As a result, 



their difference B 



dip 



B , , is also almost independent 



of W. Such insensitive behaviors to W are found also at 
moderate Q 2 . On the other hand, the i-slope difference 
in ^3 shows stronger IF-dependence as W increasing: 
their result decreases by a factor ~ 2 from the moderate 
fixed-target energies to the HERA collider energies. 

We can now compare the result of Bjf^ with ZEUS 
data [{| by combining with the i-slope arising from the 
two-gluon form factor as in Eq. J2J). Here, we simply set 
t = for the first term of J2Jl, and then it adds a con- 
stant factor of Bn to BjF^ as Bj/^(t = 0) = 4.0 + B d jj^. 
Using this relation, we find that the contribution of the 
cc state to the i-slope is not negligible, because the size 
can reach about 10 % of the whole at Q 2 = for the 
J I ip production. Similarly for the tp production, it turns 
out that the contribution of the cc state becomes largest 
(about 5 %) in the vicinity of Q 2 ~ 10 GeV 2 . The result 
of Bj/^,(W) is shown in Fig. 01 with the data measured 



for the muon and electron decay channels of J/tjj sepa- 
rately. The data indicate the increase of the t-slope with 
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FIG. 4: The 147-dependence of the t-slope Bj/^ for the 
elastic J/ip photoproduction. The data was measured in 
two leptonic decay channels, J/tp -> /x + fi~ , e+e" i. The 
sohd line shows our result Bj/^{t = 0) = 4 + Bjf^ in 
Eq. (J2J. The dotted line represents the fit to the data: 
B V (W) = B + 4a' P • In (VF/90 GeV) with B = 4.15 GeV -2 
and ap =0.116 GeV" 2 0. 



W and this energy-dependent i-slope fits well to the data 
as a form of B V (W) = B + Aa' F ■ In (IF/90 GeV) with 
B = 4.15 GeV" 2 and a' p = 0.116 GeV" 2 0. Our t- 
slope almost independent of IF is in agreement with the 
data only in the range of IF = 100 ~ 200 GeV. The 
reason for this behavior (solid line) is easily explained in 
our formulation: By P {W) giving only the IF-dependence 
in Eq. |[2J) has the form oc AA/A 2 in a simplest 

description, where >4(„4) denotes the imaginary or real 
parts of the amplitude with either longitudinal or trans- 
verse polarizations. Then, the amplitudes has a common 
factor xG{x) in the integrands. Usually, at HERA ener- 
gies, the low- a; behavior of the gluon distribution xG(x) 
is expected to be xG(x) ~ x~ x with A ~ 0.2, which is 
obtained from the analysis of the inclusive deep-inelastic 
scattering data. Here, x ~ M 2 /IF 2 at Q 2 = and the 
Q 2 -scale of G is set to Q 2 ~ to 2 . Following this sim- 
ple parametrization for xG(x), By P is completely inde- 

2 9 

pendent of IF. In fact, x and the scale Q + k ± have 
a, fej^-dependences in our formulation. The corrections, 
however, give only a weak dependence of By P on the 
energy IF, illustrated in Fig. 0] 

Thus, the energy-dependence of the t-slope, so called 
"shrinkage", observed by the ZEUS experiment, seems 
to need more sophisticated treatments for the IF depen- 
dence. For example, we may consider two possibilities 
as the origin of the IF-dependence, which we have not 
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taken into account here: one possible effect is the Gri- 
bov diffusion, which is known as a rapid expansion in the 
transverse size of the cc-dipole state, created from a pho- 
ton fluctuation as W increases. This phenomenon leads 
to the shrinkage of a diffraction peak, which can be inter- 
preted as an increase of the interaction radius Q ; another 
effect is the increase of the cc-dipole scattering off the pe- 
ripheral pion-cloud of a nucleon with increasing W . As 
suggested in 0, at relatively low W(~ fO GeV) such as 
the fixed-target energies, the contribution of the pions to 
the gluon distribution is not important. It is because the 
pions carry a smaller fraction of the nucleon momentum, 
so that the gluons inside the pions have much smaller 
momentum fraction than x ~ 0.1, which is the typical 
value of the momentum fraction at the energies. In this 
case, the authors of |l5j predict that the corresponding 
two-gluon form factor (J7J is close to the axial form factor 
of the nucleon with m\ ~ 1 GeV 2 . At higher collider 
energies of HERA (typically, W ~ 100 GeV), where it 
reaches x <C 0.01, one should take into account the con- 
tribution from the pion-cloud. At high W, therefore, the 
form factor is expected to approach the electromagnetic 
one with m\ g ~ 0.71 GeV 2 fl5| . 

In order to check how the mass scale changes with W in 
our approach, we try to reproduce the HERA data of the 
elastic J/tp photoproduction by changing a value of the 
mass scale. We neglect a contribution from the Gribov 
diffusion, because the pQCD analysis is inapplicable for 
a large size configuration due to the diffusion. 

Resulting ^-dependence of the J /if) cross section is 
shown in Fig. [SJ It is compared with ZEUS data Q 
for three representative ranges of W, i.e., 30 < W < 50 
GeV, 70 < W < 90 GeV and 150 < W < 170 GeV 
for the decay of J/tp, and 20 < W < 50 GeV, 

70 < W < 90 GeV and 125 < W < 170 GeV for 
the e + e _ decay. They are calculated over the kinematic 
range \t\ < 0.8 GeV -2 , in which our formulation could 
be applicable. Here, we set a parameter R = 1.2 for the 
skewedness effect of the off-diagonal gluon distribution, 
which changes only the normalization of the cross section 
and thus does not affect the dependence of cross section 
on t. Further, to increase the overall magnitude of the 
cross section, we carry out a rescaling of Q 2 in the gluon 
distribution and the strong coupling constant |l9j, fol- 
lowing the work of Q- It is also stressed that our naive 
treatment for the rescaling is almost insensitive to the t- 
slope, but it enhances only the normalization of the cross 
section about 30 ~ 40 %. 

The solid, dotted, and dashed curves represent our cal- 
culations corresponding to each W range, 30 < W < 50 
GeV, 70 < W < 90 GeV and 150 < W < 170 GeV. They 
show an exponential decrease very close to the data with 
increasing \t\, when we take mf g = 0.9,0.7,0.6 GeV 2 in 
respective W bins. Hence, to obtain a good agreement 
with the data, we find that the mass scale should decrease 
with increasing W at HERA energies. This is qualita- 
tively consistent with the suggestion of 15] mentioned 
above. 
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FIG. 5: The differential cross section of elastic J/ip photopro- 
duction, 7 + N — > J/tp + N, as a function of \t\, compared 
with ZEUS data measured in two leptonic decay chan- 
nels, J/ip — » n~ , e + e~ . The data are plotted for three 
representative ranges of W for respective decay channels, as 
written in the figure. The solid, dotted and dashed curve 
show our results calculated in three corresponding W ranges, 
30 < W < 50 GeV, 70 < W < 90 GeV and 150 <W< 170 
GeV. Then good fits to the data require that the squared 
mass scale has the values 0.9, 0.7, 0.6 GeV 2 in respective W 
regions. 



IV. SUMMARY AND DISCUSSION 

To summarize, in the LLA of pQCD we have formu- 
lated ^-dependences of the differential cross section for 
the diffractive ("elastic") charmonia (J/tp and tp') photo- 
and electroproductions at low \t\(< 1 GeV -2 ). To deal 
with the production of the radially excited 25 state (tp') 
together with that of J / tp, we have employed the charmo- 
nium LCWF, properly including the sub- leading effects 
due to Fermi motion, based on a technique developed 
in P, 0- Following QCD factorization, the two-gluon 
form factor of the nucleon is assumed to be process- 
independent and we define it as a function of only t scaled 
by the relevant squared mass in a dipole form. By assum- 
ing an exponential form of the differential cross section, 
following the standard experimental fit so far, we have 
calculated the diffractive i-slope, By(t), over the \t\ range 
of less than 1 GeV 2 . 

The dipole part, By V , shows that the Q 2 dependence 
presents a distinct difference between J/tp and tp' , reflect- 
ing the properties of wave functions, i.e., a node effect 
for the 25 state. The former monotonically decreases 
with Q 2 , while the latter has a negative value at small 



Q 2 and rapidly increases with Q 2 , keeping Bjf^ 



> B 



dip 

4>> 



in the realistic Q region of HERA. Then, the difference 
Bm — B^i shows a strong Q 2 dependence that tends to 
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almost zero with increasing Q 2 associated with squeezing 
of the dipole, after achieving a maximum of 0.53 GeV -2 
at Q 2 = and IV = 85 GeV. This value at Q 2 = 
is consistent with HERA data. It is noted that, to a 
good approximation, the difference should be free from a 
contribution of the Gribov diffusion, because such a con- 
tribution should be incorporated into the initial photon 
part in QCD factorization and it is expected to contribute 
at the same level for J/ip and ip'. 

The TV-dependences of By p indicate almost a flat be- 
havior at Q 2 — (or also moderate Q 2 ), contrary to a 
shrinkage of t-slope observed at ZEUS. At small x, mainly 
two effects are expected to be the origin for such In- 
dependences of the t-slope: one comes from the dipole 
part, i.e., the Gribov diffusion; another is due to the tar- 
get nucleon, i.e., the effect of hard probe scattering off 
the peripheral pion-cloud of a nucleon. We incorporated 
the latter phenomenon through the change of the mass 
scale appearing in the nucleon form factor, by fitting to 
ZEUS data for differential cross section of J /ip photopro- 
ductions at several TV ranges, although we neglected the 
former phenomenon for simplicity. To get a good fit to 
the data, we find significant decrease of the mass scale 
with increasing IV, e.g., 0.6 < m| g < 0.9 GeV 2 in the 
region of 20 < TV < 170 GeV. This mass scale decreas- 
ing with TV will provide us with information concerning 
gluon distributions of the nucleon through a probe of cc- 
dipole. Further detailed study on it might be useful for 
investigating not only the gluon distribution of the pion- 
cloud on the surface of nucleon, but also the non-linear 
property of gluons inside the nucleon, which will become 
more important with higher TV. 

On the other hand, the difference of i-slopes shows 
a behavior almost independent of TV. This quantity 
is most likely free from the Gribov diffusion, as men- 
tioned above. Unfortunately, the TV-dependence of 
has been not yet observed so far. However, its obser- 
vation has some importance and we would request early 
realization of such an observation. If the experiments ob- 



serve somewhat the TV dependence for the difference of 
i-slopes, contrary to our prediction, then we might guess 
the following two origins: 

1. The Gribov diffusion might introduce a large size 
effect upon the cc-dipole at high TV, thereby influ- 
encing the final hadronization dependence. It will 
give important information about the mechanism 
of the diffusion. 

2. Otherwise, we could doubt the hypothesis for a uni- 
versal form of nucleon form factor based on the fac- 
torization, which is very challenging. 

Finally, we stress that the finite size effect of the cc- 
dipole is not negligible for the i-slope of charmonia elec- 
troproductions at the HERA energies and Q 2 < 30 GeV 2 . 
Combining both data from J/ip and ip' productions, we 
can investigate the pure contribution of the dipole part in 
detail. Such a precise study will finally make it important 
to extract the detailed information of the gluon distribu- 
tion in the nucleon by making use of the diffractive char- 
monium photoproductions. Recently, a fascinating study 
on the energy-dependent i-slope was reported in |20j , us- 
ing the dipole saturation model with an impact parame- 
ter dependence. A part of the TV dependence could also 
arise from the non-linear dynamics of the saturation for 
the gluon distribution in the target nucleon. The analy- 
sis including such an effect we consider to be interesting 
and is planned for one of our future works. 
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